Statistical-mechanical formulation of Lyapunov exponents 
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We show how the Lyapunov exponents of a dynamic system can in general be expressed in terms 
of the free energy of a (non-Hermitian) quantum many-body problem. This puts their study as 
a problem of statistical mechanics, whose intuitive concepts and techniques of approximation can 
hence be borrowed. 



I. INTRODUCTION 

Lyapunov exponents are an important tool for the 
characterization of dynamical systems. Their very def- 
inition has a strong statistical mechanical flavor, as it 
involves 'extensivity in time' of certain quantities regard- 
less of the initial ('border') conditions. It is then natural 
to seek to express them in an explicitly statistical me- 
chanical way, in terms of a partition function. 

In this paper, building upon previous work, mainly 
of Graham Q and Gozzi p|, we show how the study 
of Lyapunov exponents, as well as Ruelle's angle, can be 
cast as a quantum many-body problem — in fact a rather 
standard one, except that the Hamiltonian is in general 
non-hermitian. As we shall see, this does not bring a 
miraculous solution to all calculational problems, but is 
serves two purposes: 

• Because the problem is formulated as a standard 
quantum many-body one, all the tools developed in 
that wider context are available. Some, but not all 
of them have already been used as approximation 
schemes for Lyapunov exponents, but others such 
as Hartree-Fock, dynamic mean-field theory, and 
the renormalization group seem promising. 

• More important, general theoretical results are re- 
expressed in a language that is often familiar. A 
typical example is when one asks whether finite- 
dimensional systems have a limit Lyapunov density 
function p(A): In this setting the question becomes 
whether a quantum finite-dimensional system has 
an extensive free energy for all chemical potentials. 
Although this does not in itself prove the existence 
of p(A) in the thermodynamic limit, it renders it 
intuitive and acceptable, at least up to the level of 
rigor of theoretical physics. 



Consider a general dynamical system: 



/i(x,?y) 



1. 



and, in particular, the version with additive noise: 
&i = /i( x ) + rji 



(1) 



(2) 



where rji is a Gaussian white variable with variance 2T^. 
In the limit of zero noise we have a standard dynamical 



system and for a particular form of / a Hamiltonian sys- 
tem. Clearly, an alternative way to study (0) is to go to 
the Fokker-Planck (or 'Kramers', or 'Liouville', depend- 
ing on the context) description of evolution of probabili- 
ties: 



P(x,i) = -H FP P(x,t) 



where 



H 



FP 



d 

dxi 



T^- + ft 

OX; 



(3) 



(4) 



Here and in what follows summation of repeated indices 
is assumed, unless otherwise stated. The Fokker-Planck 
operator acts on the space of functions of the coordinates 
x: it resembles a Schrodinger operator, although it is in 
general non-Hermitian. The noiseless limit is subtle, and 
is the subject of ergodic theory. 

Introduce now two sets of fermion a\ , b\ and boson a\ , 
1, N); and the correspond- 



fi\ creation operators (i 
ing vacuum, defined by: 



i \-)=0 i \-) = O ; Vz 



(5) 



It will turn out that all the information we search is ob- 
tained directly from the following generalization of the 
Fokker-Planck operator: 



He 



H F p + V k i{-x)(a\ai 



where we denote: 



Vfcj(x) 



dfk 
dxi 



■a\ai 



PlPi) (6) 



(7) 



Hs acts on the product space of functions of the coordi- 
nates and of the number of fcrmions and bosons of each 
type. (Clearly, Hs coincides with Hpp when restricted 
to the zero fermion and boson subspace.) This is the 
many-body system mentioned above. 

If the Xi are lattice variables, and the fi are short-range 
interactions, then the system (^|) defines a quantum (non 
hermitian) theory, also having short-range interactions. 
If the system is instead off-lattice, the Xi describe the 
position of the particles and the (ai, 6j, CKj, /3j) play the 
role of 'spin' degrees of freedom carried by the quantum 
particles: if the fi are short-ranged, then both the direct 
and the spin-spin interaction are also short range. 
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An alternative strategy, that leads to the generalized 
Lyapunov exponents, is based on replicas of fermions a] , 
b]^ with 7 = 1, q and the operator: 

H q = H FP + V kl {*){atcq + bl%1) (8) 

The Lyapunov exponents are defined as follows. The 
separation of two infinitesimally close trajectories (evolv- 
ing under the same noise, see Arnold ||) x(t) and 
x(t) + y(t) is described by the evolution of vectors in 
the tangent space: 



Vk{t) = C/fei(x o ,t)2/i(0) 



(9) 



where x Q is the initial condition and U is defined as the 
solution of the linear equation: 



U ki {t) = V kj {x,t)U ji {t) 



(10) 



The N Lyapunov exponents A^ measure the rate of 
growth of the volume in the tangent space. We construct 



A(t) = U{t)U\t) 



(11) 



and consider its eigenvalues A\{t) > A2 (t)... > AN(t) (a 
set per initial condition and/or noise realization). The 
Lyapunov exponents are 



A 4 = lim (log A t 

t — >oo t 



(12) 



Throughout this paper averages (•) are over the noise 
and/or the initial condition. Note that the Xi are auto- 
matically labeled in a decreasing order Ai > A2 > ... > 
A at. Their existence has been proved under very general 
conditions by Oseledec 

One can in fact consider the generalized Lyapunov ex- 
ponents A| (GLE), of interest as a measure of the inter- 
mittency (see Benzi et al. ||) defined for example as: 



for systems with one degree of freedom. Later, Gozzi 
jjj derived a generalization to several degrees of free- 
dom. His approach differs from ours in two respects: 
first, his expression, if averaged over initial conditions or 
noise, yields the annealed quantities. Second, and also 
important, he considers the eigenvalues of U rather than 
Ulfl , a different (though potentially interesting) quan- 
tity. In the present treatment we derive an expression 
for the quenched exponents directly by introducing aux- 
iliary fermions and bosons: the technique of supersymme- 
try H 0])- We also give expressions for the generalized 
exponents for even q. In all cases we consider the opera- 
tor UU\ rather than U: this is the reason why we need 
to introduce the fermions and bosons in pairs (a*, fo) and 
(cfr. il). 

Lyapunov exponents appear naturally in the context 
of localization problems for quantum disordered systems: 
one considers exponents associated with the product of 
transfer matrices, which play the role of the evolution 
operators, (see e.g. j|, ^|) Thus, the Green- function for- 
malism has been treated with supersymmetric (SUSY) 
techniques to extract correlations (see Balents and Fisher 
JlOf , where there is also a detailed analysis of the super- 
symmetry group and representation, and Markos |ll| for 
a path-integral formalism). In these works supersym- 
metry was used in order to obtain the Green function, 
and then the quenched quantities were obtained via the 
replica trick. 

This paper is organized as follows: In section || we 
present the formalism for the expression of usual and 
generalized Lyapunov exponents. For Hamiltonian sys- 
tems we show that the Lyapunov pairing rule follows triv- 
ially from a particle-hole symmetry, and we also discuss 
the computation of Ruelle's angle. In section III we dis- 



cuss the statistical mechanics of Lyapunov exponents of 
macroscopic systems, in particular the existence of a Lya- 
punov density function. We then present two examples of 
application: in [V to problems of Random Matrices, and 
in section^ to a Hamiltonian Mean Field model (HMF). 



X 2q 



lim - log (.4? 



(13) 



(a more standard definition will be given below in (p_7[)). 

Adopting the nomenclature of disordered systems, we 
shall call the true Lyapunov exponents 'quenched', and 
the q = 1,2,... set the 'annealed' exponents. The 
quenched exponents can be formally obtained as the an- 
alytic continuation q — > (the replica trick), but we shall 
also consider a direct evaluation here. 

Years ago, Graham 0] constructed a representation of 
the annealed Lyapunov exponents using fermions valid 



II. THE FORMALISM 



A. Basic quantities 



The average expansion of p-dimensional volumes evolv- 
ing with (||) in the space of phases can be expressed, in 
terms of U as: 



J 



det 



(14) 
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where £,h--4i P are an sets of p orthonormal vectors. In- 
troducing the cumulative Lyapunov exponents Aj as 



where T denotes time-order and 



one has 



A, = Ai + ... + A, 



A p = lim —(hiRp(x ,t)) 

i— >oo Zt 



(15) 



(16) 



The expression ( p"?j| ) is not suited to be written as an 
integral of an exponential (a partition function) because 
of the logarithm. This is the usual quenched versus an- 
nealed problem in disordered systems: we are interested 
in the average of the logarithm and not the logarithm of 
the average. Below we shall overcome it by means of su- 
persymmetry, but we shall first also give expressions for 
the annealed quantities. 

The generalized Lyapunov exponent (GLE) (see (^]) 
can be used to measure intermittency (i.e. rare tra- 
jectories having unusual Lyapunov exponents) and are 
relevant quantities in the so-called "thermodynamic for- 
malism" of chaotic systems. They are defined via the 
moments of the R p (x Ql t): 



A 2q = \\ q 



A? 



lim - \n{{R p { Xo ,ty)) 

t—tCO t 



(17) 



To the extent that the moments (R p (x ,t) q ) are suffi- 
cient to reconstruct the distribution law of R p (x 0l t) and 
the averages for all real q, one can use them to find the 
quenched average ( |l6| ) . In such cases the quenched quan- 
tities can be extracted from (|lj]) using the replica trick: 



A, 



(18) 



B. Generalized exponents 

Let us first obtain an expression for the annealed expo- 
nents \f q for integer q. This calculation is the closest to 
the construction of Gozzi |l[] , with the important differ- 
ence that we consider the limit eigenvalues of the matrix 
A = UW and not of U itself. 

Introducing p pairs of fermions ai and bi one can write: 



Rp (Xq , t) 



V4l 



(19) 



Hrit) = Vi^M^+blb,) (20) 
tp p = — a il ...a lp b ll ...b ip (21) 



The explicit time-dependence of Hi , is given by the evo- 
lution of x Eq. (||) via Vy (x). More generally, one can 
write, in terms of q replicas of the set of fermions a] and 
b': 



R q ( Xo ,t) = {-\r p T 



X Hl{x{t'))dt> 



n 



ft i 



(22) 



where 



Hf(t) = VWv(a?a} + b?b'r) 



li 7=l 



E 

i\,...,i p 



(23) 
(24) 



(see appendix ^). This is not the final expression, since 
we have not yet imposed that x evolves according to the 
equation of motion. 

The use of two sets (a's and b's) of fermions are needed 
in order to follow the rate of growth of the norm of vec- 
tors in the tangent space and not the evolution of the 
vectors themselves: we are interested in the eigenvalues 
of U'U and not in those of U. (This is obviously not a 
problem in dimension one Q ) . In dimensions larger than 
one, Gozzi Q] introduced one family of fermions and thus 
studied the eigenvalues of U : these do not in general co- 
incide with the Lyapunov exponents (see references 
and |ll| for a discussion), although they surely give rel- 
evant information. In section [V we will explicitly show 
in an example how the a and 6-fermions interfere in a 
non-trivial way. 

We now write an expression for the average over noise 
and/or initial conditions of ( p2] ) using the information 
that x(x Q , t) evolves according to (||), and the probability 
density follows (^J) . One way to do this is to express the 
weight of each trajectory as: 



(^(xo,i)) 



Dz<y(a;-x(x ,77,t / )) HV$ 



f* H?(x(t'))dt' 



e Jo 



9t| 



(25) 



where Dx means the flat functional integral over trajectories and the delta-function imposes that x(t) satisfies the 
equation of motion. This passage is just the standard textbook exercise of going from a Langevin to a Fokker-Planck 
description (see for example ) . The result is that the probability evolves through the Fokker-Planck equation, and 
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we have: 

(jg(xo,t)) = (1| ® (-|^ [e-* ff «] <|-> ® |x > (26) 



where ff g = Hpp — Hi is given in (Q). The ket (1| is the 
flat measure (l|x) = 1 and (x|x ) = ^(x — x ) is the dis- 
tribution associated to x Q . Note that the time-ordering 
is automatic, as the evolution of x is taken care of by 
Hpp- In the limit t — > oo the logarithm of (i?*(x ,i)) is 
given by the smallest eigenvalue of H q in the subspace not 
orthogonal to the vectors (1| ® (— and ^\— ) <8> |x Q ). 

We can now describe a practical algorithm for com- 
puting the GLE's. First we have to identify the smallest 
invariant subspace containing the states appearing to the 
right and to the left of (26). Clearly, the fermion num- 



bers N2 = £ z 



,7t„7 



and = J2i are conserved 



and the calculation of (i?^(x Q , t)) involves working in the 
subspace N£ — NJ = p. Furthermore, because H q and 
commute with the operators V 1 defined by: 



n 



ValV^ = -6? 



(27) 



we should look into the subspace of eigenfunctions sym- 
metric under exchange of a's and 6's, having eigenvalue 
one under the V 1 : here is where two families of fermions 
get mixed. Because the quantity (|26| ) is by construc- 
tion positive, one can see that the eigenvalue having the 
smallest real part within this subspace has zero imagi- 
nary part. 



1. Variational calculations 

A variational approach can immediately be imple- 
mented for the calculation of the generalized (and, with 
an act of faith, the usual) Lyapunov exponents. Since 
one is looking for the lowest eigenvalue of H q within a 
given subspace, one can use a variational trial function 



\<P) 



~ minRe 



{<f\<p) 



(28) 



If the family of variational functions is parametrized for 
every i, q, one obtains an approximation which is an ex- 
plicit function of q: one can then envisage computing also 
an approximation for the usual Lyapunov exponents via 



C. Ordinary (quenched) Lyapunov exponents 

In this section we construct an expression for the 
(quenched) Lyapunov exponents using a supcrsymmctry 
formalism. For the purposes of clarity, we shall do so 
in two steps: first a naive calculation that does not take 
into account the convergence of the sums, but is closer 



to the standard supcrsymmctry treatments in other con- 
texts. Next, we use the Borel transform technique to 
work more properly. 



1. Careless calculation 

We introduce a set of fermion and boson operators as 
in (^|) and (0), and define the number operators: 



N, 



bos 



N a = a\a k 
N a = a\a k 
N a + N fj 



N b = b\b k 



N a + N b 



N= N fer +N b OS 



(29) 



which will commute with all other operators. We shall 
also need: 



f = a\b\ ; f = 4pl 



(30) 



Let us now introduce, for a given trajectory x(t), the 
quantity: 

fff (t) = Vw(x)(<4ai + b\h + a\ai + fifi) (31) 
With these notations we compute Z(/j,, p) as: 



(32) 

an one can easily show that (Appendix H) : 



det[l 



A] 



det[l 



-4] 



(33) 



where A(t) is associated with the trajectory x. This func- 
tion will generate all the Lyapunov exponents for the tra- 
jectory as: 



G(ji) 



lim 

f^oo t 



dp 



N 



lim > — 

t-voo 1 



^ 1 + e-^Aj 

J — l J 

(34) 

since for large t, G(fi) is a ladder with steps at the val- 
ues of /I corresponding to the Aj's (i.e. it is the integral 
of the Lyapunov distribution function). Just as in the 
previous subsection, we wish to calculate (G(/i)). Again, 
this is directly done in the Fokker-Planck formalism by 
the quantity: 



\ W t-^oo t dp 



(35) 



fi=fj, 
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with (Z[fx, (J,)) written as: 

( ( f >L \e- t ^ Hs+ ' lN f^ + ' iNb ^(-l) NB '"\(l) R ) (36) 

where Hs = Hpp — Hf is given in (^). The left and right 
eigenvectors are 



\4>r) 



dX„ P(x r 



J + f\ 



(37) 



where P(x Q ) is the initial condition distribution. We can 
also write (G(p)) as an expectation value: 

(G( M )) = (Me- t{Hs+ ^ } N fer (-l) N »°*\<t> R ) (38) 

Equation (B3[) can be expressed as a formal series in 
powers of exp[—t[j,) and exp(— tp,). This same series is 
reproduced by (p^), or, in averaged version, (^). Each 
term of the form exp[— t(nifi + n2p)] corresponds to an 
expectation value in the subspace of n\ fermions and n 2 
bosons. Clearly, ([33]) has a very small convergence radius: 
this is because the number of bosons, unlike the number 
of fermions, is unlimited. In other words, expressions 
(|33|), (||^) and ( |36"1) have only a formal meaning. One can 
still work with them if at the end of a derivation one can 
resum exactly the series, in which case one has in fact 
performed the analytic continuation in /i and /2. 



2. Borel transform 

Let us now give a more proper construction to the su- 
persymmetric expression of the Lyapunov exponents. It 
will turn out that the formalism that emerges is not much 
more complicated than the one in the previous para- 
graphs. 

Briefly, the Borel transform technique consists of going 
from the formal series: 



h(y) — a Q + aiy + a 2 y 2 + ■■■ 
to its convergent transform 



(39) 



«i 



"2 ,2 



h (y) =a + —y+—y 



I!" 2! 

which can then be inverted. In our case, we shall take 



(40) 



h(y) ~ a + aiy + a 2 y 2 + 



det[l + ye'^A] 
det[l + ye-^A] 



(41) 



(cfr. Eq. (|33j)), and define as the Borel-transformed par- 
tition function: 



Z B (^) = h B (y)\ =1 



(42) 



that is, we are dividing by (ni + n 2 )\ the term having 
ni fermions and n 2 bosons. Repeating the construction 
above, it is easy to see that: 



{z B {^n)) = {<% | e -*(**+^+^)(-i)^|<^) 



(43) 



which is just like ( pq ) except that the left and right eigen- 
vectors are 



dx Q P(x )5(/ + /)|-)®|x ) 
(-|®(-|5(/ t + / t ) (44) 



where the function g(x) is: 

oo 

9(x) = E 



n!v n' 



We can now define as before: 



(G B (n)) 



lim 



ld(Z B (^ft)) 



dfi 



(45) 



(46) 



which we can also write as an expectation value: 

(G B ( M )) = (<fi | e -*{H s +„iV} iV /er (-l)^|^> (47) 

Interestingly enough, we can retrieve the information di- 
rectly (without the need to antitransform) from G B (/i), 
since one can show (Appendix |^) that: 



JV 



3=1 

again gives the same ladder function. 



(48) 



3. Symmetries 

Although we shall not make here much use of the sym- 
metries in the problem, let us briefly mention them for 
completeness. Clearly, the Hamiltonian Hs is invari- 
ant under all the (supersymmctric) transformations ro- 
tating simultaneously the cti , f3i , ai , hi so as to leave the 
quadratic form / + / and the fermion-boson part of Hs 
invariant. Then, expectation values can be written in the 
standard form: 



(O) =Tr[(-l) N f° r CO] 



(49) 



with C supersymmetric: 

C = e -^s + (M+-)(^ + ^}| B )(0 B| (50) 

Both the original and the Borel transformed versions 
have the same symmetries. 

One can show that supersymmetry is responsible for 
the fact that Z(/j,,fj,) and Z b (/j,,ij,) are independent of 
/i: the constancy of the normalization is indeed the un- 
derlying reason why we can use the method to obtain 
quenched averages. (See Ref. 0) for a discussion in 
detail of the supersymmetry group and representations) . 
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D. Discussion 



the Vij read, 



The calculation of generalized exponents \f q is done, as 
we have seen, by computing the lowest eigenvalue within 
a subspace of the Hilbert space. This is because the 
large-i limit automatically projects onto the correspond- 
ing eigenstate. Approximate and numerical methods for 
the estimation of the ground state of Schrodinger-like op- 
erators abound in the literature, we have already men- 
tioned the variational principle. 

If one wishes to extend these results for the quenched 
exponents, an analytic continuation to q — ► is needed. 
This is easily done (although with a leap of faith) when 
an explicit expression for all even q is available. Such will 
be the case in a variational, a perturbative or a mean-field 
computation. 

On the other hand, the supersymmetry method yields 
the Lyapunov exponents without the need of any con- 
tinuation. However, there is a price to pay: expressions 
( |38| ) and ([l?1 ) involve a sum of terms within subspaces 
of any number of bosons. In the Borel transformed ver- 
sion (|47| ) this sum is convergent for all finite t. However, 
one can see that the largest term corresponds to a bo- 
son number of the order of exp(Ait), and this number 
grows as we consider larger times. In other words, we can 
only perform this sum for finite t, and only then make 
t — > oo. Again, this is no problem if an analytic expres- 
sion is available (perhaps as a result of an approxima- 
tion), but it does seem problematic to attack a problem 
numerically this way. 



E. Hamiltonian systems: Pairing rule and Ruelle 
Angle 

Damped Hamiltonian systems are a particular case of 
the dynamics (|2|), which can be written as 

dH 

opt 

P = « 7Pi + 7 2 rii ( 51 ) 

oqi 

where 7 measures the intensity of the coupling to the 
bath. The simplectic structure of Hamilton's equations 
has consequences. 



1. Pairing of exponents 

One of the proprieties we can easily infer from our 
formulation is the pairing- rule for the Lyapunov systems. 
This pairing rule was proved for a class of dynamical 
systems (named 'quasi- hamiltonian' by Dressier fl5|| ); in 
our formalism it is the result of a particle-hole symmetry. 
This is most easily seen for the annealed exponents. In 
the case 



V 



H 



2 





d 2 V 
dqidqj 



(53) 



Denoting x% — qi for (i = l,...,N) and Xi = pi for 
(i = N + 1, 2N), we consider the transformation: 



where 



a i = e H a 3 >>, 



-1 

1 



(54) 



(55) 



(here and in what follows and 1 are the null matrix 
and the identity matrix in the N x N space). 

Under this transformation the quantities from jl9| ) 
transform as: 



H q ^ H g -2qj{N-N fer ) 
Vv|0) = f/>N-p\0) 



(56) 
(57) 



This particle-hole symmetry, a consequence of the sym- 
plectic structure of the evolution operator, was noted by 
Gozzi [0. One then gets: 

A^=A^_ 9 + {N-p)2 r , (58) 

which in terms of the exponents becomes 

\%> = -2 qi - Aj* (59) 



Using (|18| ) and ( |15|) we can infer that the Lyapunov ex- 
ponents are symmetric around — 2 , Note that this result 
holds for arbitrary r\i (t) , and not only in the case in which 
it is a white noise. 



2. Ruelle's angle 

The difficulty in computing the Lyapunov exponents 
is due to the fact that, in the tangent space, any instan- 
taneous frame generically turns. In order to define an 
angle associated with these transformations Ruelle 
defined a rotation number (or Ruelle frequency -RF- af- 
ter Dressier |]l7[ ) and proved an additive theorem for it. 
To construct this number we can use a polar decomposi- 
tion of the evolution operator: U — A? Q with A = UU^ 
simplectic and symmetric and Q simplectic and orthog- 
onal: 



Q 



X Y 

-Y X 



(60) 



V(q) 



(52) 



Clearly, T = X + iY is unitary; and we can extract a 
rotation number from det(T) — e" n (the determinant of 
an unitary matrix is a pure phase). 
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The Ruelle angle can be extracted directly from U (see 
fyl) in the following way [^9|: write the evolution oper- 
ator as 



U = 



U A U b 
U c U D 



(61) 



Then, m(t) will be the argument of z(t) — detlUA+iUs)- 
We can write: 



2im(t) _ 



z{tf 



det(U A + iU B f 



z{t)z{t)* det(U A U A + U B U B ) 



(62) 



where we have used the fact that, since U is a real sim- 
plectic matrix, U A U B = U B U A . 

We can use fermions to write the determinant quanti- 
ties as: 



N 



N 



det(U A +iU B ) = (-\l[a qi Te tv n^+^JH ( 63 ) 



i=i 



(here we adopt the notation a qi a Pi for the fermion cor- 
responding to the i — th position or momentum degree 



of freedom (i = 1,...,N). We can further use the b 
fermions to obtain the square of this determinant. To 
compute the denominator of (|62|), we first note that 
O = U A U\ + U B U B is the top-left block in the A = UW . 
To have his determinant in the denominator we can use 
the relation: 



(_| e « 9s /5 9s Te tV (-l) Nh °'e a *iP«* +a »t p *i\-} 

(64) 

But, since O is an N x N block of A, it will generically 
have a projection in the direction of the N largest eigen- 
values, so for large times: 



det(l + O) 



det(l + 0(t)) ~ det(0(t)) ~ e 



s > t 



(65) 



where 5*> is the sum of the largest N Lyapunov expo- 
nents. Putting everything together, we can write 



N N 

(e Mm ®) = (l|®<-|e^ n°* 6 « Ii« +i4 i )(^ ! +4 I )(- 1 ) JVboae< ^ +< ^|-)®|xo) (66) 

k=l 1=1 



Note that the quantity we have calculated is not quite 
the quenched averaged angle, but rather the average of 
the exponential. 



III. MACROSCOPIC SYSTEMS 

An active field of research is the information that Lya- 
punov exponents can provide in extensive systems in 
thermodynamic limit. The natural object to study is 
the Lyapunov density function: 



N 



or, better, the cumulative version: 



C(n) 



p(\)d\ 



(67) 



(68) 



In particular C(0) is the sum of positive exponents, re- 
lated to the Kolmogorov-Sinai (KS) entropy. The exis- 
tence of a thermodynamic limit for the Lyapunov den- 
sities has been conjectured (|lj| ^0|). A problem imme- 
diately arises: the Lyapunov exponents are themselves 
the result of the limit t — > oo, and the question as to 



whether this, and the thermodynamic limits commute is 
not obvious. Typical examples when they do not is when 
there are macroscopic motions that take times that di- 
verge with the size. We shall find a clear example in the 
Hamiltonian mean- field model below (Section |y|): the 
particles organize in an anisotropic object which can turn 
collectively like a rotor: In contact with a bath there is 
diffusion of the collective angle, but, since the moment of 
inertia scales with N, the collective motion is absent if we 
consider N — > oo before t — > oo. Sinai |2l[| has shown that 
in this last order of limits the densities are well defined in 
a system of confined particles with pairwise interactions: 
we shall see that this result is very natural. 

In the previous sections we have shown that C(/i) is 
the large time limit of Ct(/i) with 



Ct(fi) = G s ( M ,t) 



1 d\nZ B (p,fl,t) 



dji 



(69) 



(the logarithm in the r.h.s. has no effect, since the 
normalization is one). 

Consider this expression: Z B is a partition function 
associated with the 'quantum' hamiltonian Hs, where 
the time plays the role of an inverse temperature, /x of a 
chemical potential and Ct(fJ-) of the derivative of a free 
energy density with respect to the chemical potential (i.e. 



8 



a particle number per unit volume). The Lyapunov den- 
sity is hence a form of compressibility. 

As mentioned in the introduction, if the original prob- 
lem is on a lattice and has nearest neighbor interactions, 
the fermions and bosons are also lattice variables inter- 
acting with the nearest neighbor variables (through Vij ) . 
We have then a 'quantum' lattice problem with short 
range interactions. On the other hand, the system could 
be off-lattice, and the Xi be a set of d-dimensional vec- 
tors describing the position of the particles interacting 
via short-range pair forces /i(x) = — Xj). The 

variables (ai,bi,ati, Pi) play the role of 'spin' degrees of 
freedom carried by quantum particles, both the direct 
and the spin-spin interaction are also short range. 

All in all, we are asking whether a quantum theory with 
short range interactions has a good thermodynamic limit 
with a well defined free-energy density. There is only one 
non-standard feature if we ask for the t — > oo taken before 
the thermodynamic limit: this is, as we have seen, like 
asking in a statistical mechanical problem about the zero 
temperature limit taken before the thermodynamic limit 
- sometimes a tricky question. 

The arguments on extensivity become more subtle if 
we wish to study the thermodynamic limit of the largest 
Lyapunov exponent: this is like asking in a particle sys- 
tem not what is the chemical potential needed to create 
a certain particle density, but rather to create a single 
particle in the whole system: clearly this is a question of 
order 0(l/N). We shall return to this point in section 

An interesting special case is the behavior of the Lya- 
punov exponents close to zero in a system with soft 
modes. In the present context this concerns the prop- 
erties of a statistical mechanic problem around /i = 0, 
i.e. free from external chemical potential. For example, 
from Eq. (^8|) we have: 

(p(A - 0)> = {<pf \e-* Hs N fer (N fer + N bos )\^) (70) 

A very intriguing possibility that immediately comes 
to mind when working in the present framework is that 
of studying universality properties in critical points using 
renormalization-group ideas and techniques. 



IV. RANDOM MATRICES 

In this section we use our formalism to derive some 
results already obtained for some Random Matrix (RM) 
models. For brevity we shall only do this in the pure 
fermion (annealed) case, although the supersymmetric 
approach can also be applied. 

The main lesson we shall obtain is that these systems 
become, by virtue of the disorder, interacting fermion 
problems. As such, they can be very well attacked by 
some of the many methods devised for such cases: Feyn- 
man diagrams of course, but also ressumations such as 
Hartree-Fock. 



A. Weak disorder expansion 

The first model where we can show the power of our 
approach is the one proposed by Derrida & al.||. They 
study the weak disorder expansion of the quenched Lya- 
punov exponents for a product of the form: 

P = II U * U t = B + eB (71) 

t 

where B a is a fixed matrix, B is a random matrix, and e 
is a small parameter. We shall study the case with B a 
multi-dimensional, Gaussian white noise, with zero mean 
(as the finite mean can be safely included in the constant 
matrix) : 

(fly(i)) = 

(B^it^Buih)} = <7 l3M 5(h-t 2 ) (72) 

Following Derrida we study the case with the matrix B 
having non-degenerate eigenvalues, well separated. 

The first step in order to use our formalism is the set- 
up a continuous-time variant of the problem: 

U(dt) = 1 + B dt + eBdt (73) 

Next, we exponentiate this expression. Due to the non- 
continuous character of the random term B the correct 
form of the exponential is: 

U{dt) = e Badt+Bdt-$(B*)dt* (?4) 

This means that the evolution of U will be given by (as 
in Eq. ©): 

(t) = B olJ + By (i) - i ( B « B ^ > ( 75 ) 
k 

We can derive now the path-integral form of R\ as: 
<-| J] ...67 |e-«*| fl a]l..a]X...b]l\-) 

7—1 7=1 

(76) 

(7 is the replica index). We used as a Hamiltonian 

H = V ij (x)J2(aV*] + b?b]) (77) 

7=1 

The average (Rf) will be expressed by integrating the 
gaussian noise in J76| ) which will transform H into 

H' = H + e 2 H! (78) 
H = B oll j2{aVa< + bVb1) (79) 

7=1 

7', 7=1 

+bVb^bZb] + a^b^bla]) (80) 
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We have used the fact that we can diagonalize B a and 
work on that base (which will not change the fermionic 
states). We have also arranged the creators and the de- 
structors of fermions in normal form. The expression (|80|) 
can be proved to be the correct one by going back to the 
Suzuki- Trotter product, performing the averages there, 
and then reconstructing the continuous version. We are 
now in the possession of a time-independent Hamilto- 
nian. Note that after integration of the noise we get 
a non-trivial result only because of the presence of two 
types of fermions: without them we would have lost the 
higher moments of the noise. The replica trick is easy to 
implement for perturbation expansions, since these yield 
an explicit dependence as a polynomial in q of each term 
of the expansion. This is exactly the program we carry 
on now. First, the eigenstates of H Q are: 



also Crisanti p3[ ) in order to mimic some systems that 
show strong chaos. A continuous-time version of this 
model, in the case of only one spatial degree of freedom 
(N = 2) is a linearized evolution of the form 



V(t) 



1 

v(t) o 



(85) 



where 77 is an Gaussian noise with a mean r and a devi- 
ation a. 

Using this definition we can compute the average of 



(JJ«(t)) = / P(r,)d V (r\T 



e \l H{ V {t'))dt' 



|^t) (86) 



where 



(81) 



and the corresponding eigenvalue is: 



(82) 



(ej are the eigenvalues of B Q ). The first non-zero order in 
e is in e 2 ; in order to compute it we must use first-order 
perturbation theory to obtain: 

WoWMt) = - { -^M(2q) + 3q 2^1lM. (83 ) 

For the second order perturbation theory (which will 
give the contribution in e 2 ) we must identify the states 
connected by the perturbation. We can now retain the 
terms linear in q, and hence obtain the quen ched average 
(the coefficients of the linear term (see (|17|)): 

E 2 {BijBji) 



P P P 

- 4 EEEE 

i=l 1=1 j = l m>p 
4 P P 



(BijB mi ) (BjiBi, n ) 



Si - e. 



\ ^ \ ^ \ ^ \ " (Bi m B j n ) {B m ^B n j) 
i=l j = l n>p m>p J 

This expansion is the continuous-time, Gaussian white 
noise equivalent of the expression obtained by Derrida|p|. 



B. Parisi-Vulpiani 



(87) 



7=1 



and 



9 



// : : EK^P + b V b p) + »?(*) EKW + b V h V 



7=1 



7=1 



This Hamiltonian form of the exponential quantities as- 
sociated with the Lyapunov exponents can be checked 
to be true using small-time developments in the Suzuki- 
Trotter formula, just as in the previous section. We can 
integrate the noise in d8(]) and obtain: 



(R(t) q ) = (ip q \T e tH 



(89) 



where the averaged Hamiltonian can be again checked in 
the discretized version: 



H ' = EW +a p+ 6 « t ^)+ r "EK t <+ 6 P t ^) 

7—1 7—1 

2 



+ y(EK t < + ^) 



(90) 



Let us first concentrate on the zero average case r = 0. 
We make a redefinition of the fermions as: 



,t 



a„ — > ii a 



i« 9 



dp — > t2dp : a 
1 



t 2 



t ^ > 

(91) 



Another model that can be revisited is the random 
matrix model introduced by Parisi and Vulpiani p2f (sec 



In terms of the new fermions, the dependence on a be- 
comes explicit: 
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H' 



+ b1%) + \ \Y^a\ + b^ q ) 

7=1 \7=1 / 



(92) 



We now have only to diagonalize this time-independent 
Hamiltonian on subspaces with exactly two fermions of 
each replica. 

We have done this numerically for values of q up to 25 
within the 3 9 -dimensional basis generated by: 

a?b$ ; a?bf ; ^=«&? f + a?b?) (93) 



(symmetric under Eq. (p7|)). For small q it is easy to 
diagonalize H, for example for q = 1 we obtain \\ = 

(2<7 2 )i ' 

Let us note here that the Hamiltonian (|9(]) is, for the 
case q = 1 the matrix 

V2\ 

a 2 \/2r (94) 
v \/2> V2 / 

essentially the same used by Anteneodo and Vallejos 
to derive the same results. 

Let us now turn to the case in which the noise has 
non-zero average. We can transform ( p0[ ) using transfor- 
mations © with 



h r 



(95) 



The Hamiltonian becomes H' = ^JrH" with 



H" = ^(^t a 7 + 6 7t 6 7 ) + ^ (a 7t a 7 + 6 7t 6 7 ) + 
7—1 7—1 

+ S 2 (^2( a P a «+ b V b V^J (96) 

where and s 2 = -§73 . There is a crossover (see Lima |24| ) 
in the Lyapunov exponent dependence between the limits 
of small and large s (see Fig. [l]). 

Finally, let us give an example on the role played by 
spatial structures. We consider a one dimensional sys- 
tem with a noise 77 (Eq. (|35"1)) random variables with 
correlations depending on the distance between sites: 

(%•(*)) = 

(Vij(t)Vki(t')) = o 2 {\i-]\){5 lk 5 3 i + 5 a 5 ]k )5{t-t') 

(97) 

We shall compute the annealed Lyapunov exponent (A 2 ). 
The Hamiltonian is 



M 



M 



1=1 



Lm=l 




— l-(s A 2)/8 

■ Numerical 

■■■ -0.31*s A (2/3) 

■ - annealed 

- - replica trick 



s A 2 



FIG. 1: The greatest Lyapunov exponent versus the strength 
of the noise (s=cr) for a symplectic random matrix (the mean 
of the noise is 1). At a >> 1 replica approximation and the 
power law fit are presented. The dashed line is the annealed 

A 2 

Lyapunov exponent (-£■). For small s , the Lyapunov expo- 

(s) 2 

nent is given by 1 — ^-g-; the inset shows the check of this 
formula against the numerical data. 



After averaging it becomes 

M 

H = Z)(4 i a f ( + 6 g,V) + l^^ 2 (0)(4 i a 9! +6| ! 6 9 J 2 



M 



1 = 1 



M 



^ a (\l k\) ( ^ + h ^ + a ^ Qgi + ^tjs 

k^l=l 



(99) 



The system becomes explicitly translational-invariant. 
We have to diagonalize this Hamiltonian on the states 



1 M 

— al, b\. t 



1=1 



M 



(100) 



(98) 



where M is the total number of sites. Within this sub- 
space, H is identical to the matrix (p3) with r = and 
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we get: 



a 2 = AiC£<j 2 (*) + 2<7 2 m = Mj^imm) (ioi) 



i=l 



again a result obtained before in Refs. [E2l 



We can use now the function defined in (|32|) in order to 
compute the Lyapunov spectrum. The supersymmetry 
Hamiltonian becomes a sum of single particle terms: 



V. HAMILTONIAN MEAN FIELD SYSTEM 

The Hamiltonian mean field system we consider p5[ is 
composed of N coupled rotators with a classical hamil- 
tonian: 



2 N 



^[1 - cosfe - qj)] 



(102) 



This system has (in the canonical ensemble) a phase tran- 
sition at T c — 0.5 J between a ferromagnetic and a para- 
magnetic phase |25|, |2(| |27j 

The equations of motion are, if we allow for noise and 
dissipation: 



dqi 
dt 
dpt 
dt 



-JM X shift + J My cos ft - 7p l + \f 7 2rffr\ i 



(103) 



where we have introduced: 



M x 



TV ^ 



cos ; M y 



losing, (104) 



H 



(i) 



—— [JM X sin q l - J M y cos ft 



dpf 



+jv l p \ v l v% + J(M X cos ft + My sin qi)v l p \ v\. 

-Jv p \ cos q,C l - Jt^ sin q,S l (107) 
(here no summation over i is implied). For compactness 
we have defined v l for I — 1, 2, 3, 4: 



and the collective operators: 



(108) 



The matrix Vij governing the evolution operator in 
tangent space is: 



with 



V = 



1 

-E -7I 



(105) 



COS qi 



N 



2^ sin ft (109) 



Eij = - — cos(ft - qj) for i ^ j 



N 



JM X cos qi + J My sin Qi ~~ (106) 



Not surprisingly, the operator Hs is a mean-field (quan- 
tum) operator itself, and we can solve the problem with 
any of the usual methods. For example, introducing ex- 
plicitly the collective variables using the functional delta- 
functions: 
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D[M X ] S(NM x -J2cos qi ) = 

i 

D[M y ] 5(NM y -J2 shl< l*) = 

i 

D[C 1 } 5(NC l cos *) = 

i 

D[S 1 } 8(NS l -Y,vl sin g 4 ) - 



D[M X ]D[M X ] e -M x (NM x -Y, i cosq i 

D[M y ]D[M y ] e ~ !il v< NM v-T.i«™ii) 
D[C l ]D[C l ] e- d '' wc '-^< cos ^ 
D[S l ]D[S l ] e -^ NS, -^*< sin «*) 



(110) 



we can write: 



D[M X ]D[M X ] H D[C l ]D[C l ]D[S l ]D[S l ] e~ N ^ dt (M^+M y M y+ c'c l +s l s')~w) 
i 

I 



(111) 



where W is the action for a single pair of variables (q, , pi 



e w = (^ITe-'W^^+^^-l)^!^) 

(H2) 

(no summation). The left and right eigenvectors are 



l^) = / da i P(a; i )e' W,+ ^ 5t |->®|z<> 



(113) 



where for simplicity we have assumed that all pairs of 
variables (qi,Pi) have the same initial distribution, and 
hence left and right vectors are in product form. The 
single particle effective Hamiltonian (H^^) equals: 



H 



M x cos qi — My sin qi —Cv cos qt — S l v l q . sin qi 

^ ' ' ' (H4) 

which has in principle a time dependence through the 
collective variables. We wish to evaluate < Z(fi, fj,') > 
by saddle-point method in the thermodynamical limit. 
The saddle point equations for the ordinary collective 
variables read: 

9W w dW 
M x = ■ 7 ^- M T = 



dM x 
dW 
dM v 



dM x 
dW 



(115) 



We also have saddle point equations for the sixteen 
fermionic and bosonic collective variables C l ,C l ,S ' ,<S\ 



We have 20 equations of type (115) which we can solve 
assuming that the system is at t = already in thermo- 
dynamic equilibrium - the vector to the right in ( ]112| ) is 
Gibbs-distributed. This implies that the collective vari- 
ables may be constant in time, an assumption we verify 
later. First of all, it is easy to see that the saddle-point 
values of the C l ,C l ,S l and S L vanish (Appendix |E|) . Next , 



one can see that if the endpoint of the trajectories is left 
free (i.e. we are not conditioning to a specific arrival 
point), then causality implies that M x = M y = (Ap- 
pendix |e|). Under these assumptions, the only variables 
with non-zero saddle-point values are: 



M x 



{cos qi) M 



(sin q l ) M 



(116) 



where the average (•) m is taken with the single particle 
dynamics: 



dq 

lit 
dp 



-JM sin g — 7P + y2rfTrj 



(117) 



Indeed, this yields the equilibrium value of M, as solution 
of the equation: 



M = 



h(y) 
h{y) 



is 



MJ 



(118) 



(where / are the Bessel functions, see [p5|). What we 
have shown is that the Lyapunov density function of the 
system is, to leading order in N the sum of Lyapunov 
densities for a single particle system moving according to 
Eqn. 0j\ ). 

In the limit of zero coupling to the bath 7 = the 
exponents for a single particle are zero, as a consequence 
of conservation of energy and the pairing rule. Note that 
the fact that the energy of each particle is conserved sep- 
arately is an artifact of the large iV-limit: there is a cou- 
pling between the fluctuations at the following order in 
1/N. In conclusion, in the thermodynamical limit the 
function G(/u) is a step of height ~ 27V at zero. 

A vanishing largest Lyapunov exponent Xi has already 
been obtained both numerically ([g(| |27j) and analyti- 
cally ([^3| |28[ ]) in the paramagnetic phase, but not in 
the ferromagnetic phase. This is not in contradiction 
with our results: according to our calculation one can 
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1,5- 



0,5- 



■■■N=100 
■-N=500 
-N=1000 




-0,03 -0,02 -0,01 0,01 0,02 0,03 
x 



FIG. 2: The microcanonical Lyapunov spectra of the HMF 
system with N particles at the energy E — 0.3 presented as 
^ffl (see equation $5%) for different numbers of particles (N). 
Only the first half of the spectra is computed, the other half 
being obtained by the pairing rule. 
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FIG. 3: L yapu nov exponent for one particle with Kramers 
dynamics (117) versus the energy for 7 = .01 (dots) and 7 = 
.001 (squares) respectively. E = 0.75 is the energy at the 
transition point (for J = 1) and the Lyapunov exponent is 
automatically zero above this point. The exponents are scaled 
with 73 . 



VI. CONCLUSIONS 



still have a vanishing fraction of non-zero Lyapunov ex- 
ponents, that do not contribute to the density function 
in the large- N limit. If we wish to calculate the largest 
exponent Ai we should take our calculation to the next 
order. 

Figure || shows a numerical calculation of the spectrum 
of Lyapunov exponents for various values of the number 
of particles N. The simulations are carried out in the 
microcanonical ensemble at an energy where the greatest 
Lyapunov exponent does not show an important varia- 
tion with N (E = 0.3, see Even though the com- 
puter time requirements do not allow us to get close to 
the thermodynamical limit, these numerical results show 
a reasonably good agreement with the hypothesis of a 
step spectrum. 

At finite 7 the Lyapunov distribution function is no 
longer a step, even in the thermodynamic limit. We 
can obtain it very easily by solving the one-particle sys- 
tem (117) numerically. The results depend strongly of 7; 
there is a scaling law in 7 1 / 3 below the critical tempera- 
ture and an identical zero Ai in the paramagnetic phase 
(see Fig. |). 

Surprisingly, after rescaling the largest Lyapunov expo- 
nent behaves very much like the one computed by Firpo 
(|p8|) in the canonical ensemble. There are also qualita- 
tive and quantitative resemblances with the one obtained 
in micro-canonical simulations of a large (but obviously 



Writing the Lyapunov exponents as a statistical me- 
chanic object provides a different perspective of the prob- 
lem, allowing to transpose much of the knowledge and 
intuition developed in that wider context. 

Perhaps the most clear example are the questions re- 
lated to the thermodynamic limit. We have seen that 
the existence of a limit Lyapunov density function is the 
kind of extensivity property that most theoretical physi- 
cists would accept without proof — at least for the times 
that do not diverge with the system size. 

As to the different approximation schemes, a case in 
point is the continuous product of random matrices Sec- 
tion IV . We have seen there that the problem is mapped 



finite) number of particles (see 
elude that, at finite 7, G(fi) shou 



2f| |27jj). We can con- 
d have two sigmoids of 



height ~ N located symmetrically around zero. 



into a system of interacting fermions. We have treated it 
as a perturbative expansion of weak interaction reobtain- 
ing the results of Ref. ||. However, once in the language 
of interacting particles, one immediately thinks of other, 
more global approximations; we have already mentioned 
Hartree-Fock, which can easily be implemented both in 
the pure fermion or in the supersymmetric formalism. 
Similarly, the standard mean-field treatment of the model 
of section [v] can be extended for the calculation of the 
exponents. This also suggests that it may be interesting 
in general to construct a local mean-field approximation 
for problems with space: this could give a simple ana- 
lytic handle on the spatial structures involved with each 
exponent. 

A question into which we have not looked in detail 
is intermittency. The generating function of generalized 
Lyapunov exponents X i q involves a distribution function 
in space that depends on q. Physically, this arises be- 
cause we are conditioning the probability of a trajectory 
to having an unusual value of the exponents. Hence, 
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studying the lowest eigenvector of H q gives us informa- 
tion on the spatial structures responsible for intermittent 
behavior. 

Yet another interesting question is to look at systems 
at or near criticality, and borrow methods and ideas from 
the rich theory of critical phenomena to infer results on 
the behavior of the exponents there. 
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APPENDIX A 

In this appendix we show how the use of fermions 
transform the quantity R defined in ( p"4| ) can be writ- 
ten as a certain matrix element of the hamiltonian (|2~T 
(equation [l9|)). 



We put: 



Let a\ be a set of fermion operators 
[a,;, a]] + = % 
We can encode ( |T4| ) by writing: 

0{t)=a\y k {t) 



(Al) 



(A2) 



6(t) = [Vualai ,0(t)] 



(A3) 



It is easy to check that the yt then satisfy (Q). The 
solution to ( |A3| ) is: 

0(t) = (Te tVkia ^ O(0) (Te tVk ''' a l' a ' , y 1 (A4) 
which implies: 

f Te tv kia {ai\ t (Te tv "'" a l' a "y 1 ^U Jt a] (A5) 

Let us now add a second set of fermions b k , and con- 
struct the operator: 



so we can write UjiU r iOjbl as: 



(A6) 



(note that the same constructions can be made with 
bosons instead of fermions.) Denoting |— ) the fermion 
vacuum, Eq. (A7) implies: 



{-\ah-a>i b h ...b ip \T 



^Vij{a\aj+b\bj) 



U klh ...U kv i v U jlh ...U 3v i v (~\a il ...a lp al i ...al v \-) 

I 



(A8) 



Summing over ii, i p the first part transforms into: 

U klh U Jlh ...U kplp U Jplp = (UU% in ...(UU^ kpJp (A9) 

Denoting with j p ) (fci, ...,k p )) the signa- 

ture of the permutation betwe en t he two sets of indices 
the second (fermionic) part of (A8) equals 



a{{i u ...,i p ) ->■ (ki, -,k p ))a((ji, ...,j p ) -> (i u ...,i p )) 

(A10) 

Using ( |A9| ) (|A8| ) and the properties of the signature it 
follows that (A8) equals 



(pO A jiki-Aj n h p v((ji,-,jp)^(ki,...,k p )) 

Jljfel 5 jpi^p 

(All) 



(remember that A = U'U). After summing on fci, k p 
the expression ( All[ ) is, by definition: 



(p!) 2 det \(yl, ,yt )t [tfu] (y ii: ,y ip )] (A12) 



APPENDIX B 



In this appendix we deri ve f ormula (|32j ) . Denoting 
the fermion vacuum, Eq. (A7) implies: 



(Bl) 



and: 



(_| e / + (T e tVkl ( a k ai+h * h ^-» N t°rt\ e /|_) = 

(_| e ^«* e -^/ar e 7A,>atbt|_^ _ det [ 1 + e -Mt^] (B2) 
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The last equality is most easily checked by performing a 
rotation of the fermions to diagonalize A. Then we have: 

(_| e 6ia 4e -/iJV / „* e AiraJit|_j = 



(_| e b*< e -^ e ^a; t 6; t |_ ) = 

n fc (_| e fcM-e- f 'V^V'|-) 



(B3) 



where A are the eigenvalues of A (no summation on k); in 
the last expression there is no sum on i in the exponents. 
Each factor can be developed into 1 + e~ Mt Ai, so the 
whole expression will be: 



JJ(1 + e-^Ai) = det[l + e'^A] 



(B4) 



Consider now the same steps, but now replacing the 
fermions a\ and b\ by bosons a\ and (3\, as in (po|); then 
we obtain the analogue of (B2): 



(-|e /f ^re tVkl{a ^ ai+ ^ l}l) ~ tlNhost ^ e*\-) 



= det[l-e~ tlt A]- 1 (B5) 



which can again be proved by performing a rotation o f the 
bosons to diagonalize A. After the same step as in ([B3|) 



the terms that survive are those with an equal number of 
the creation and destruction operators and we can use 



(-|(a*/3t«)|->=n! 



to obtain the analogue of (B4) 



-lit A, 



-2 fit j^L 



II: 



- e-^Ai 



(B6) 



(B7) 



One can change the sign in (B7) we can introduce some 
"negative norm" bosonic states, or as: 



1 



det\l + A] 



(_| e / _t (Ye n/fci(a fc ai+/ ^ /3i) ) e~f\-) 



(-\e? fr^'t^^+flft)) (_i)JV»« e /|_) (B8) 



APPENDIX C 

In this appendix we show the ladder structure of the 
functions G defined in (^) and (|47]) . The common ingre- 
dient after developing the exponentials and matching the 
terms allowed by the conservation of bosons and fermions 



J 



(_|(/t -|- yt)" ^7-(_l)^ O a e tV^^4 a l+ 6 * 6 0+*^*i( Q; t Q; (+/ 3 lA)-Mt(^!>os+Ar/er)^ jV/ er (/ + /)"|-) (CI) 

We can commute the "evolution" operator with Nf er and apply it on the right ket; the action of the operators number 
of bosons and fermions Nb os + Nf er (on the left ket) can be easily computed so the expression becomes: 

e - n » tH (P-Pr NF (M4 + M0tr H (C2) 

The expansion of the powers in terms of individual creators an destructors: 

«-W £ (-in "*» N >" n^) tj ^ l-> (03) 



Here there is no sum on repeated index inside the parenthesis. Expanding further Nf er in terms of fermions operators 
this term can be written as: 

e-^*(n!) 2 V V ( I {hm - Piai)kl a ' a ' + b ' bl A k > ( °' °' + TT ^ ~ A* ^ + I TCI) 

But for each I, at j ^ I and fcj > a factor of the product above will be: 

{ _ Aj)kl f utfcpapfipn ^(-i^-^r^r^ 1 ^" 1 !-)^^^^^!-) ^ (C5) 

and we can easily see that this equals zero. So, inside each term of the sum on 1, only that with ki = n and the 
others kj — w ill survive. 
The sum (C3) simplifies to: 

^(_l)»-1^2/_|^-l n-l tn-l^tn-li_w_| 6 oWtfe Q t 6 t|_\ 
e- n ^(n!) 2 J2— " W UlU. = e-n^Y^Afi-l) 11 - 1 (C6) 

I n - I 



16 



Now we can reconstruct the sum on n for the different 
cases; the Borel construction gives uses the series 

N 00 r_p-M* a \ n N 



j=l n=l 



3=1 



The series from ( C7 ) are convergent at each time and for 
each trajectory (in a well behaving, smoot h p otential); 
one can formally sum the general terms (C6) without 
the —7 and obtain: 



N 



N 



j=l n=l 



3=1 



(C8) 



Using the exponential form of Aj, in the limit t — ► 00, 
both (C7) and (|C8|) will give 



N 

g(m) = g b (m) = E^ a -^ ( C9 ) 

3=1 

APPENDIX D 

In this appendix we derive the fourth term in e of the 
weak disorder expansion treated in section ( I V A ) . All 
the eigenstates of H Q can be obtained by destroying and 
creating fermions in ip a . Due to the fact that the per- 
turbation contains two creators and two destructors the 
states connected by the perturbation with i\> are only 
those that differ from it by one or two fermions. The 
"replica composition" of the states connected by the per- 
turbation must be the same (two different replicas are not 
connected). 

One set of states is formed by destroying one fermion 
from the first p and creating one (from the last N — p); 
both fermions must be in the same replica family: 

j < p m > p (Dl) 
These eigenstates will contribute to the Lyapunov ex- 



ponent by a second order perturbation term: 

(lptrfjm\V\ipo){^o\V\ipb^jm) 



EEE 

7=1 j=l m>p 
q p 

EEE 

7 — 1 j—l m>p 



Sj S r, 



|v|Vo>M>M^ 



(D2) 
(D3) 



The linear part in q of this term is: 

(BjjB m j) (BjiBim) 



p p p 



(D4) 



i— 1 l — l j — l m>p J 

The states that have two different fermions are of three 



types: 



(D5) 



i < P j < P m > p n > p 
and each of them will contribute with terms of the form 

(ip\vw )(Mvm 



EEE EEE 



Ej + Si 



(D6) 



7—1 j—l m>p j' — l i—1 n>p 

The linear term in q will be: 



P P 



<?EEEE- 

i—1 j — l n>p m>p 



+ 6i 



(D7) 



APPENDIX E 

In this appendix we study the consistency equations 
for the collective variables in the HMF model. As we 
can see in from the equation ([34]) we need quantities at 
fi = £l; the derivation implied by the definition of the G 
function will be carried out on < Z(fj,,ji) > after the sad- 
dle point evaluation. So, we must take the derivative of 
the expression (111). In this point, the derivative reads: 





1 






- ~ ~Z~s 





r 



E 



dX k 



dfi 



(El) 



where we take Xk, for k = 1, 20, to be the vector of the 
collective v ariab les, and Z s is the value of < Z(n, fi) > at 
the saddle (111). We see that the second part of the rhs 
of this equation is zero because of the consistency equa- 



tions ( |115| ); in conclusion we are interested in the values 
of the collective variables at /1 = /2 and only the direct 
dependence on the two variables of W will be interesting 
for us. With this assumption we write the consistency 
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equation for the boson and fermion variables as: 



S l = 





Te 




sin qi(- 










e w 








(4 ] 


Te 


Pi 


singi(- 










e w 








($\ 


Te 




cos g^- 


i)^:i 








e w 








_ (#1 


Te 


Pi 


cosgi(- 


i)<l 





(E2) 



(E3) 



(E4) 



(E5) 



All those equations contain expectation values of single 
boson and fermion operators which necessarily vanish. 
This is normal, as we could expect from the beginning 
that those variables integrate out, contributing only with 
a non-exponential prefactor. This prefactor (as expected 
from the supersymmetry considerations ) will be one in 
the limit fi = p, but in the process of derivation will give 
terms of order O(j^). 



There are four equations left can be divided in two 
parts. The first parts contains the variables M x and M y : 



M r = J- 



(E6) 



M x = -J- 



g| T e-^+^J, sing, - sin g ^z4)(-l)^%W) 



(E7) 



The first term of the rhs of these equations is zero (as demanded by the causality) and the second is also zero, due to 
the mismatch in the fermion-boson operators. 

Finally, the last two equations are the only nontrivial ones: 



^ ) |Te- t ^// + ^cos( gz )(-l)^|^ ) ) 



M _ (0W|re- t TO + ^Bin( g Q(-l)<^W) 



Now the hamiltonian H^ff is much simpler: 



d S 2 d 



(E8) 
(E9) 

(E10) 



r 



where we used also the rotational symmetry (which define where GW (^t) characterize the Lya puno v spectrum of a 

a saddle manifold) on the space of M x and M y which single particle with the dynamics ( 117 ); this spectrum 

allow us t o fix My — and M x = M; looking back at contains in fact only two exponents. 
(Ill) and (El) we can conclude that 



n \ at r 1 dw ^v) 
G(m) = N hm — w 



t— >oo te 



(Ell) 



Formula ( |112| ) combined with the simple expression of 
i?e// (E10) allow us to infer that 

G(fi) = NG^ip) (E12) 
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